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Abstract 

Schubert polynomials were discovered by A. Lascoux and M. Schiitzenberger in the study of 
cohomology rings of flag manifolds in 1980’s. These polynomials generalize Schur polynomials, 
and form a linear basis of multivariate polynomials. In 2003, Lenart and Sottile introduced skew 
Schubert polynomials, which generalize skew Schur polynomials, and expand in the Schubert 
basis with the generalized Littlewood-Richardson coefficients. 

In this paper we initiate the study of these two families of polynomials from the perspective 
of computational complexity theory. We first observe that skew Schubert polynomials, and 
therefore Schubert polynomials, are in (when evaluating on non-negative integral inputs) 
and VNP. 

Our main result is a deterministic algorithm that computes the expansion of a polynomial / 
of degree d in Zpi,..., Xn] in the basis of Schubert polynomials, assuming an oracle computing 
Schubert polynomials. This algorithm runs in time polynomial in n, d, and the bit size of the 
expansion. This generalizes, and derandomizes, the sparse interpolation algorithm of symmetric 
polynomials in the Schur basis by Barvinok and Fomin (Advances in Applied Mathematics, 
I8(3):271”285). In fact, our interpolation algorithm is general enough to accommodate any 
linear basis satisfying certain natural properties. 

Applications of the above results include a new algorithm that computes the generalized 
Littlewood-Richardson coefficients. 


1 Introduction 

Polynomial interpolation problem. The classical polynomial interpolation problem starts 
with a set of data points, (ai,6i), ..., (a„+i, bn+i), where ai,bi G Q, Oj / aj for i ^ j, and asks 
for a univariate polynomial / G Q[x] s.t. /(oj) = for i = 1,... , n -|- 1- While such a polynomial 
of degree < n exists and is unique, depending on different choices of linear bases, several formulas, 
under the name of Newton, Lagrange, and Vandermonde, have become classical. 

In theoretical computer science (TCS), the following problem also bears the name interpola¬ 
tion of polynomials, studied from 1970’s: given a black-box access to a multivariate polynomial 
/ G F[xi,... ,Xn] (F a field), compute / by writing out its sum-of-monomials expression. Several 
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algorithms have been proposed to solve this problem [Zip79, BOT88, KY88, Zip90, KSOl]. A nat¬ 
ural generalization is to consider expressing / using the more powerful arithmetic circuits. In this 
more general setting, the problem is called the reconstruction problem for arithmetic circuits [SYIO, 
Chap. 5], and the sum-of-monomials expression of / is viewed as a subclass of arithmetic circuits, 
namely depth-2 SII circuits. Reconstruction problems for various models gained quite momentum 
recently [BC98, Shp09, SV08, KS09, AMSIO, GKLll, Kayl2, GKL12, GKQ14]. 

As mentioned, for interpolation of univariate polynomials, different formulas depend on different 
choices of linear bases. On the other hand, for interpolation (or more precisely, reconstruction) of 
multivariate polynomials in the TGS setting, the algorithms depend on the computation models 
crucially. In the latter context, to our best knowledge, only the linear basis of monomials (viewed 
as depth-2 Eli circuits) has been considered. 

Schubert polynomials. In this paper, we consider the interpolation of multivariate polynomials 
in the TGS setting, but in another linear basis of multivariate polynomials, namely the Schubert 
polynomials. This provides another natural direction for generalizing the multivariate polynomial 
interpolation problem. Furthermore, as will be explained below, such an interpolation algorithm 
can be used to compute certain quantities in geometry that are of great interest, yet not well- 
understood. 

Schubert polynomials were discovered by Lascoux and Schiitzenberger [LS82] in the study of 
cohomology rings of flag manifolds in 1980’s. See Definition 8 or Definition 12 for the definition^, 
and [Mac91, ManOl] for detailed results. For now we only point out that (1) Schubert polynomials 
in Z[xi,... ,Xn] are indexed by v = (ui,... ,u„) G N", denoted by Yv;^ (2) Yv is homogeneous of 
degree YJi=i 

Schubert polynomials have many distinguished properties. They form a linear basis of multi¬ 
variate polynomials, and yield a generalization of Newton interpolation formula to the multivariate 
case [LasOS, Sec. 9.6]. Also, Schur polynomials are special Schubert polynomials ( 13 (1)). A Schu¬ 
bert polynomial can contain exponentially many monomials: the complete homogeneous symmetric 
polynomials are special Schubert polynomials ( 13 (2)). It is not clear to us whether Schubert poly¬ 
nomials have polynomial-size arithmetic circuits. Because of these reasons, interpolation in the 
Schubert basis could not be covered by the aforementioned results for the reconstruction problems, 
unless the arithmetic circuit complexity of Schubert polynomials is understood better: at present, 
we are only able to put Schubert polynomials in VNP. 

While Schubert polynomials are mostly studied due to their deep geometric meanings (see e.g. 
[KM05]), they do have certain algorithmic aspects that have been studied shortly after their intro¬ 
duction in 1982. Indeed, an early paper on Schubert polynomials by Lascoux and Schiitzenberger 
was concerned about using them to compute the Littlewood-Richardson coefficients [LS85]. That 
procedure has been implemented in the program system Symmetrica [KKL92], which includes a set 
of routines to work with Schubert polynomials. On the other hand, the complexity-theoretic study 
of the algorithmic aspects of Schubert polynomials seems lacking, and we hope that this paper 

^Definition 12 is one of the classical definitions of Schubert polynomials, while Definition 8 defines Schubert 
polynomials in the context of skew Schubert polynomials. 

^In the literature, it is more common that Schubert polynomials indexed by permutations instead of N". These 
two index sets are equivalent, through the correspondence between permutations and N" as described in Section 2. 
We adopt N" in the introduction because they are easier to work with when dealing with a fixed number of variables. 

^Algorithms in this work run in time polynomial in the degree of the polynomial. By (2) this is equivalent to that 
the indices of Schubert polynomials are given in unary. 
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serves as a modest step towards this direction. 


Our results. Our main result is about deterministic interpolation of sparse polynomials with 
integer coefficients in the Schubert basis, modulo an oracle that computes Schubert polynomials. 
The complexity is measured by the bit size of the representation. 

Theorem 1. Suppose we are given (1) black-box access to some polynomial f G Z[xi,... , Xn], 
/ = ^.^gpOvTv, Ov / 0 G Z with the promise that deg(/) < d and |r| < m; (2) an oracle that 
computes the evaluation of Schubert polynomials on nonnegative integral points. Then there exists 
a deterministic algorithm, that outputs the expansion of f in the basis of Schubert polynomials. The 
algorithm runs in time polynomial in n, d, m, and log(X^ygp kvD- 

In fact, Theorem 1 relies on the algorithm in Theorem 9 which applies to a more general 
setting: it only requires the linear basis where the leading monomials are “easy to isolate.” See 
Our techniques for a more detailed discussion. 

Theorem 1 generalizes and derandomizes a result by Barvinok and Fomin, who in [BF97] present 
a randomized algorithm that interpolates sparse symmetric polynomials in the Schur basis. As 
mentioned, Schur polynomials are special Schubert polynomials, and the Jacobi-Trudi formulas 
for Schur polynomials yield efficient algorithms to compute them. So to recover Barvinok and 
Fomin’s result, apply our interpolation algorithm to symmetric polynomials, and replace the ifT 
oracle computing Schubert polynomials by the efficient algorithm computing Schur polynomials. 
Likewise, for those Schubert polynomials with efficient evaluation procedures^, we can get rid of 
the ffV oracle to obtain a polynomial-time interpolation algorithm. 

Our second result concerns the evaluation of Schubert polynomials. In fact, we shall work 
with a generalization of Schubert polynomials, namely skew Schubert polynomials as defined by 
Lenart and Sottile [LS03]. We will describe the definition in Section 2. For now, we only remark 
that skew Schubert polynomials generalize Schubert polynomials in a way analogous to how skew 
Schur polynomials generalize Schur polynomials. A skew Schubert polynomial, denoted by Tw/v! is 
indexed by v < w G N™ where < denotes the Bruhat order^. Schubert polynomials can be defined 
by setting w to correspond to the permutation maximal in the Bruhat order. 

Theorem 2. Given v, w G N™' in unary, and a G N” in binary, computing is in #P. 

Corollary 3. Given v G N"" in unary, and a G N"" in binary, computing Tv(a) is in #P. 

Note that in Theorem 2 we have v,w G N™', while in Corollary 3 we have v G N”. This is 
because, while for Schubert polynomials we know W for v G N” depends on n variables, such a 
relation is not clear for skew Schubert polynomials. 

Finally, we also study these polynomials in the framework of algebraic complexity. 

Theorem 4. Skew Schubert polynomials, and therefore Schubert polynomials, are in VNP. 

^For example, there are determinantal formulas [ManOl, Sec. 2.6] for Schubert polynomials indexed by 2143- 
avoiding permutations {$i < j < k < £ s.t. a{j) < a{i) < a{i) < cr(fc)), and 321-avoiding permutations {$i < 
j < k s.t. cr(fc) < a{j} < fy{i)). 2143-avoiding permutations are also known as vexillary permutations and form a 
generalization of Grassmannian permutations. 

^Bruhat order on codes is inherited from the Bruhat order on permutations through the correspondence between 
codes and permutations as in Section 2. 
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Applications of our algorithms. A long-standing open problem about Schubert polynomials 
is to give a combinatorially positive rule, or in other words, a algorithm for the generalized 
Littlewood-Richardson (LR) coefficients, defined as the coefficients of the expansion of products of 
two Schubert polynomials in the Schubert basis. They are also the coefficients of the expansion 
of skew Schubert polynomials in the Schubert basis [LS03]. These numbers are of great interest 
in algebraic geometry, since they are the intersection numbers of Schubert varieties in the flag 
manifold. See [ABS14, MPP14] for recent developments on this. 

The original LR coefficients are a special case when replacing Schubert with Schur in the above 
definition. It is known that the original LR coefficients are ^^P-complete, by the celebrated LR rule, 
and a result of [NarOG]. Therefore the generalized LR coefficients are also #P-hard to compute, 
while as mentioned, putting this problem in ^^P is considered to be very difficult ~ in fact, we were 
not aware of any non-trivial complexity-theoretic upper bound. On the other hand, by interpolating 
skew Schubert polynomials in the Schubert basis, we have the following. 

Corollary 5. Given w, v E N™', let P = {u G N™' | 7 ^ 0}. Then there exists a deterministic 

algorithm that, given access to an #P oracle, computes (aw”|u G N™) in time polynomial in |w|, |r|, 
and log(X;ugp “w'")- 

The algorithm in Corollary 5 has the benefit of running in time polynomial in the bit size of 
Therefore, when |r| is small compared to our algorithm is expected to lead to a 

notable saving, compared to those algorithms that are solely based on positive rules, e.g. [KogOl]. 
([KogOl] furthermore only deals with the case of Schubert times Schur.) Of course, in practice we 
need to take into account the time for evaluating Schubert polynomials. 

In addition, we note that Barvinok and Fomin’s original motivation is to compute e.g. the 
Littlewood-Richardson coefficients, Kostka numbers, and the irreducible characters of the symmet¬ 
ric group. See [BF97, Sec. 1] for the definitions and importance of these numbers. Since our 
algorithm can recover theirs (without referring to a oracle), it can be used to compute these 
quantities as well. Note that our algorithm is moreover deterministic. 

Our original motivation of this work was to better understand this approach of Barvinok and 
Fomin to compute the LR coefficients. This topic recently receives attention in complexity theory 
[NarOG, MNS12, B113], due to its connection to the geometric complexity theory (GOT) [Mulll, 
MNS12]. Though this direction of generalization does not apply to GCT directly, we believe it 
helps in a better understanding (e.g. a derandomization) of this approach of computing the LR 
coefficients. 

Our techniques. We achieve Theorem 1 by first formalizing some natural properties of a linear 
basis of (subrings of) multivariate polynomials (Section 4.1). These are helpful for the interpolation 
purpose. If a basis satisfies these properties, we call this basis interpolation-friendly, or I-friendly 
for short. Then we present a deterministic interpolation algorithm for I-friendly bases (Theorem 9). 
We then prove that the Schubert basis is interpolation-friendly® (Section 5). 

Technically, for the interpolation algorithm, we combine the structure of the Barvinok-Fomin 
algorithm with several ingredients from the elegant deterministic interpolation algorithm for sparse 

®While the proofs for these properties of Schubert polynomials are easy, and should be known by experts, we 
include complete proofs as we could not find complete proofs or explicit statements. Most properties of Schubert 
polynomials, e.g. the result. Lemma 14, Proposition 17, can also be obtained by a combinatorial tool called RC 
graphs [BB93]. We do not attempt to get optimal results (e.g. in Corollary 16 and Proposition 17), but are content 
with bounds that are good enough for our purpose. 
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polynomials in the monomial basis by Klivans and Spielman [KSOl]. We deduce the key property 
for Schubert polynomials to be I-friendly, via the transition formula of Lascoux and Schiitzenberger 
[LS85]. The concept of I-friendly bases and the corresponding interpolation algorithm may be of 
independent interest, since they may be used to apply to other bases of (subrings of) the multivariate 
polynomial ring, e.g. Grothendieck polynomials and Macdonald polynomials [LaslS]. 

We would like to emphasize a subtle point of Schubert polynomials that is cruicial for our 
algorithm: for if the monomial x“ is in Y^, then v dominates u reversely; that is, Vn > 

Vn + Vn -1 >Un + Un-ij ■ ■ ■ ■, Vn + ''' + vi > Un + ''' + ui- While by no means a difficult property 
and known to experts, it is interesting that the only reference we can hnd is a footnote in Lascoux’s 
book [LaslS, pp. 62, footnote 4], so we prove it in Lemma 14. On the other hand, in the literature 
a weaker property is often mentioned, that is v is no less than u in the reverse lexicographic order. 
However, this order turns out to not suffice for the purpose of interpolation. 

Comparison with the Barvinok-Fomin algorithm. The underlying structures of the Barvinok- 
Fomin algorithm and ours are quite similar. There are certain major differences though. 

From the mathematical side, note that Barvinok and Fomin used the dominance order of mono¬ 
mials, which corresponds to the use of upper triangular matrix in Section 4.1. On the other hand, 
we make use of the reverse dominance order, which corresponds to the use of lower traingular matrix 
in Proposition 18. It is not hard to see that the dominance order could not work for all Schubert 
polynomials. We also need to upgrade several points (e.g. the computation, and the bounds on 
coefficients) from Schur polynomials to Schubert polynomials. 

From the algorithmic side, both our algorithm and the Barvinok-Fomin algorithm reduce mul¬ 
tivariate interpolation to univariate interpolation. Here are two key differences. Firstly, Barvinok 
and Fomin relied on randomness to obtain a set of linear forms s.t. most of them achieve distinct 
values for a small set of vectors. We resort to a deterministic construction of Klivans and Spielman 
for this set, therefore derandomizing the Barvinok-Fomin algorithm. Secondly, our algorithm has 
a recursive structure as the Barvinok-Fomin algorithm. But in each recursive step, the approaches 
are different; ours is based on the method of the Klivans-Spielman algorithm. As a consequence, 
our algorithm does not need to know the bounds on the coefficients in the expansion, while the basic 
algorithm in [BF97, Sec. 4.1] does. Barvinok and Fomin avoided the dependence on this bound 
via binary search and probabilistic verification in [BF97, Sec. 4.2]. However, it seems difficult to 
derandomize this probabilistic verification procedure. 

Organization. In Section 2 we present certain preliminaries. In Section 3 we define skew Schubert 
polynomials and Schubert polynomials, and present the proof for Theorem 2, Corollary 3 and 
Theorem 4. In Section 4 we define interpolation-friendly bases, and present the interpolation 
algorithm in such bases. In Section 5 we prove that Schubert polynomials form an I-friendly basis, 
therefore proving Theorem 1. We remind the reader that skew Schubert polynomials are only 
studied in Section 3. 

2 Preliminaries 

Notations. For n £ N, [n] := {!,..., n}. Let x = (xi,..., Xn) be a tuple of n variables. When 
no ambiguity, x may represent the set {xi,..., Xn}- For e = (ei,..., e„) £ the monomial with 
exponent e is := ... x®". Given / £ Z[x] and e £ Coeff(e, /) denotes the coefficient of 
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in /. X® (or e) is in / if Coeff(e,/) 7 ^ 0, and Ef := {e G N” | G /}. Given two vectors 

c = (ci,..., Cn) and e = (ei,..., Cn) in Q"", their inner product is (c, e) = Yl7=i Each c G Q"" 
defines a linear form c*, which maps e G Q”' to (c, e). 

Codes and permutations. We call v = {vi,... ,Vn) G N”' a code. We identify v with v' = 
(fi,..., Vn,0, • • •, 0) G for m > n. The weight of the code v, denoted by |v|, is ^ 

code V G N” is dominant if vi > V 2 > ■ ■ ■ > Vn- We dehne a partition to be a dominant code, and 
often represent it using a. v is anti-dominant if vi < V 2 < ■ ■ ■ < Vk and Vk+i = ■ ■ ■ = Vn = 0- 

For N G N, Sjsi is the symmetric group on [A^]. A permutation a G Sn is written as 
a{l),a{2), ... , (t{N). We identify a with a' = a{l),..., a{N),N + 1,..., M G Sm, for M > N. 
The length of a, denoted by |( 7 |, is the number of inversions of a. That is, |(t| = |{(i,i) : i < 
j-,a{i) > a{j)}\. 

Given a permutation a G Sn, we can associate a code v G by assigning Vi = \{j ■ j > 
i, a{j) < cr(i)}|. On the other hand, given a code v G N"' we associate a permutation cr G Sn{N > n) 
as follows. {N will be clear from the construction procedure.) To start, cr(l) is assigned as vi + 1. 
cj( 2 ) is the (u 2 + l)th number, skipping (t( 1 ) if necessary (i.e. if cj( 1 ) is within the first (u 2 + 1 ) 
numbers). (j{k) is then the {vk + l)th number, skipping some of cr(l),... , a{k — 1) if necessary. For 
example, it can be verified that 316245 gives the code (2, 0,3, 0,0, 0) = (2,0, 3) and vice versa. 

Given a code v its associated permutation is denoted as (v). Gonversely, the code of a permu¬ 
tation a G Sn is denoted as c(cr). It is clear that |(t| = |v|. 

Bases. Let i? be a (possibly nonproper) subring of the polynomial ring Z[x]. Suppose M is a 
basis of i? as a Z-module. M is usually indexed by some index set A, and M = {tx | A G A}. 
For example A = N”' for R = Z[x], and A = {partitions in N”} for R = {symmetric polynomials}. 
f G R can be expressed uniquely as / = ax ^ 0 G Z, tx G M, and a finite F C A. 

A construction of Klivans and Spielman. We present a construction of Klivans and Spielman 
that is the key to the derandomization here. Given positive integers m, n, and 0 < e < 1, let 
t = [m^n/e]. Let d be another positive integer, and fix a prime p larger than t and d. Now define 
a set of t vectors in as KS(m, n, e, d,p) := {c(^\ ..., c^*)} by mod p, for i G [n]. Note 

that cp^ < p = Oim^ndje). The main property we need from KS is the following. 

Lemma 6 ([KSOl, Lemma 3]). Let ...,he m distinct vectors from N"' with entries in 
{0,1,... , dj. Then 

Pr are distinct for j G [ml] > 1 — m?n/t > 1 — e 

fcek] 

On ffP and VNP. The standard definition of ffV is as follows: function / : U„gN{0,1}" ^ Z 
is in ffV, if there exists a polynomial-time Turing machine M and a polynomial p, s.t. for any 
X G {0,1}"', f{x) = \{y G {0, IjPfo) s.t. M accepts {x,y)}\. In the proof of Theorem 2 in Section 3, 
we find it handy to consider the class of Turing machines that output a nonnegative integer (instead 
of just accept or reject), and functions / : U,igN{ 0 , 1 }” ^ N s.t. there exists a polynomial-time 
Turing machine M and a polynomial p, s.t. for any x G {0,1}”, f{x) = X^^gjo i}p(") 

^This is known as the Lehmer code of a permutation; see e.g. [ManOl, Section 2.1]. 
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As described in [dCSW13], such functions are also in as we can construct a usual Turing 
machine M', which takes 3 arguments {x,y,z), and M'{x,y,z) accepts if and only ii z < M{x,y). 
Then ^ M'{x,y, z) = Y^yM{x,y). Note that z G {0, for some polynomial g as M is 

polynomial-time. 

The reader is referred to [SYIO] for basic notions like arithmetic circuits. VP denotes the class 
of polynomial families {/njnsN s.t. each /„ is a polynomial in poly(n) variables, of poly(n) degree, 
and can be computed by an arithmetic circuit of size poly(n). VNP is the class of polynomials 
{gn}nm s.t. gnixi,...,Xn) = E(ci Cm)e{0 1}^ ■ ■ ■, Xn, Cl,..., Cm) where m = poly(n), and 

{fn} nGN is in VP. Valiant’s criterion is useful to put polynomial families in VNP. 

Theorem 7 (Valiant’s criterion, [Val79]). Suppose cj) : {0, 1}* —>■ N is a function in #P/po\y. Then 
the polynomial family {/njnsN defined by fn = i}" is in VNP. 

3 Skew Schubert polynomials in and VNP 

In this section we first define skew Schubert polynomials via the labeled Bruhat order as in [LS03]. 
We also indicate how Schubert polynomials form a special case of skew Schubert polynomials. We 
then put these polynomials, and therefore Schubert polynomials, in and VNP. Also note that 
it is more convenient to work with permutations instead of codes in this section. 

Definition of skew Schubert polynomials. The Bruhat order on permutations in Sm is defined 
by its covers: for cr, vr G Sat, a < vr if (i) cr“^ 7 r is a transposition Tjfc for i < k, and (ii) |cj| -|- 1 = | 7 r|. 
Assuming (i), condition (ii) is equivalent to: 

(a) a{i) < cr{k); 

(b) for any j such that i < j < k, either a{j) > (j{k), or a{j) < cr(i). 

This is because vr gets an inversion added due to the transposition Ttk- So in no position between 
i and k can a take a value between a{i) and a{k). Else, the number of inversions in vr will change 
by more than 1. Taking the transitive closure gives the Bruhat order (<). The maximal element 
in Bruhat order is vro = V, V — 1,..., 1, whose code is d = (N — 1, V — 2,..., 1). 

The labeled Bruhat order is the key to the definition of skew Schubert polynomials. While 
naming it as an order, it is actually a directed graph with multiple labeled edges, with the vertices 
being the permutations in Sn- For a < tt s.t. a~^7r = Tst, s < j < t and b = a(s) = 7r(t), add a 

labeled direct edge as a That is, for each a < tt, there are t — s edges between them. 

For any saturated chain C in this graph, we associate a monomial where e(C') = 

(ei,... ,eAr_i), and e* counts the number of i appearing as the first coordinate of a label in C. 
A chain 

(il.^l) 02 , 62 ) Um,bm) 

(To y (Ti ^ . . . y Cm 

is inereasing if its sequence of labels is increasing in the lexicographic order on pairs of integers. 
Now we arrive at the definition of skew Schubert polynomials. 

Definition 8. Let d and ttq be as above. Given two permutations c and tt, s.t. c is no larger than 
TT in the Bruhat order (c <7r), the skew Schubert polynomial 

V/a(x) := ( 1 ) 

c 
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summing over all increasing chains in the labeled Bruhat order from a to vr. The Schubert poly¬ 
nomial := yTTo/o-- 

Now in the increasing chain (or any chain in the labeled Bruhat order), each edge increases |(t| 
by 1. So the number of edges in an increasing chain from a to vr is |7r| — |(t|. 


Skew Schubert polynomials in ffF and VNP. Before describing the algorithm for skew 
Schubert polynomials, we note the following. First, by the correspondence between codes and 
permutations, from a code v G N” we can compute (v) G Sn in time polynomial in |v|. Also, we 
have N = maxjgj^jjuj -|- i} < |vjl + n < |vo| -|- n. Second, the length of the path from cr to vr is 

Kl ~ I'^l- 

To start with let us see how Corollary 3 follows from Theorem 2. 


Proof of Corollary 3. Given v, we can compute (v) G Sm- Note that < |v| -|- n. Then form 
w = {N — 1, — 2,..., 1). Invoke Theorem 2 with (v, w,a). □ 


Proof of Theorem 2. Consider the following Turing machine M: the input to M are (1) codes 
V, w G N™ in unary; (2) a = (oi,..., a^) G N"" in binary; and (3) a sequence s of triplets of integers 
{si,ji,ti) G [A^] X [A^] X [A^] where Sj < ji < ti, N = |w| -|- m, and i G [.£], i = |w| — |v|. (Note that 
all the conditions on s can be checked efficiently.) The output of M is a nonnegative integer. 

Given the input (v,w, a, s), M uses s as a guide to compute an increasing chain from (v) to 
(w) in the labeled Bruhat order oi Sn- Let (Tq = (v). Suppose the chain to be constructed is 

ai- At step i, 0 < i < i, M also maintains an exponent vector 


o-Q 




CTl 


{hM) 


(jehe) 


e, G and the label {ji,bi). 

To start, M sets oq = (0,..., 0), and {jo, bo) = (0, 0) (lexicographically minimal). Then when 
the step i — \ finishes, M maintains e.j_i, crj_i, and (jj_i,6j_i). Then at step i, based on {si,ji,ti), 
M performs the following. First, it computes CTj = ai-iTs^ti, and checks whether |crj_i| -|-1 = |(Ti|: if 
equal, then continue; otherwise, return 0 (not a valid chain). Second, it sets e* by adding 1 to the 
jjth component of ej_i, and keeping the other components same. Third, it computes bi = (Ji{ti), 
and checks whether {ji,bi) is larger than (ji_i,6j_i) in the lexicographic order: if it is larger, then 
continue; otherwise, return 0 (not a valid chain). 

When the ilh. step finishes, M obtains cj£ and e^. It first checks whether ai = (w): if equal, 
then continue. Otherwise, return 0 (not a valid chain). M then computes as the output. 

This finishes the description of M. Clearly M runs in time polynomial in the input size. M 
terminates within I steps where I = |w| — |v|, and recall that w and v are given in unary. 

Finally note that yw/v(a) is equal to M(v, w, a, s), where s runs over all sequences of 
triplets of indices as described at the beginning of the proof. By the discussion of jfP at the end 
of Section 2, this puts the evaluation of Tw/v in n 


Proof of Theorem 4- Let us outline the proof of Theorem 4, which basically follows from the proof 
of Theorem 2. By Valiant’s criterion (Theorem 7), to put Tw/v in VNP, it suffices to show that the 
coefficient of any monomial in l^/v is in ^P. Therefore we consider the following Turing machine 
M' , which is modified from the Turing machine M in the proof of Theorem 2 as follows. Firstly, M' 
takes input (v, w, a, s), where a is thought of as an exponent vector, and is given in unary. Second, 
in the last step, M' checks whether d — equals a or not. If equal, then output 1. Otherwise 
output 0. It is clear that this gives a ^P algorithm to compute the coefficient of x®. The only 





small problem here is that in the literature ([BiirOO, Prop. 2.20], [Koi05, Thm. 2.3]) we can find, 
Valiant’s criterion is only stated for multilinear polynomials.® But it only takes a little more effort 
to overcome this; essentially, this is because the degree is assumed to be polynomially bounded, 
so a can be given in unary. First, note that N (as in the beginning of the proof of Theorem 2)is 
an upper bound on the individual degree for each Xj. Then introduce N copies for each variable 
Xi- Since a = (ai,... ,a„) is given in unary, we can assume w.l.o.g. that each a* is an V-bit string 
(ajp,... ,aj^Ar), and if a* = k, the first k bits are set to 1, and the rest 0. (These conditions are 
easy to enforce in the definition of Mb) We then use aij to control whether we have the jth copy 
oi Xi, or 1, using the formula aijXi + 1 — Uij. After this slight modihcation, the Valiant’s criterion 
applies, and the proof is concluded. □ 

4 Sparse interpolation in an interpolation-friendly basis 

4.1 Interpolation-friendly bases 

Let M = {tx I A G A} be a basis of a Z-module R C Z[xi,... ,Xn], indexed by A G A. Given a 
function K : A x N ^ M'*', M is called K-bounded, if G M, the absolute values of the coefficients 
in the monomial expansion of tx G M, tx G Z[xi,..., x„] are bounded by K{\, n).® 

For a (0,1), nonsingular matrix A, La := {c* G (Q”)* | 3c' G (Z''')"',c = Ac'}. Note that c' 
is a vector with positive integer components. Also recall that for c G Q", c* denotes the linear 
form determined by c. M is La- compatible, if for every tx G M, we can associate an exponent e;^ 
s.t. (1) for any c* G La, c* achieves the maximum uniquely at e;, over Et^ = {e G N" | G 
(2) ex / ca' for A / A'; (3) the coefficient of in tx is 1. x®^ (resp. ca) is called the leading 
monomial (resp. leading exponent) of tx w.r.t. La- By the conditions (1) and (2), the leading 
monomials are distinct across M, and for each tx, the leading monomial is unique. We assume 
that from A it is easy to compute ca, and vice versa. In fact, for Schubert polynomials, A is from 
A = Z", and e = A. 

Combining the above two definitions, we say M is (AT, Lyi)-interpolation-friendly, if (1) M is 
A'-bounded; (2) M is A^-compatible. We also call it {K, A^)-friendly for short, or I-friendly when 
K and La are understood from the context. 

• A trivial example is the monomial basis for Z[x], where K = 1, A \s the identity transforma¬ 
tion, and a leading monomial for x® is just itself; 

• For symmetric polynomials, the basis of Schur polynomials (indexed by partitions a) is 
(A', Lyi)-friendly, where (1) K{a,n) = Y^|a|! by Proposition 11, (2) A = {ri,j)ij£[n] where 
Tij = 0 if z > j, and 1 otherwise, by the fact that every exponent vector in Sq is dominated 
by a ([BF97, Sec. 2.2]). {A is the upper triangular matrix with I’s on the diagonal and above.) 
The associated leading monomial for Sq, is x". In retrospect, the fact that Schur polynomials 
form an I-friendly basis is the mathematical support of the Barvinok-Fomin algorithm [BF97]. 

®It is well-known though that Valiant’s criterion works for even non-multilinear polynomials. However, the only 
reference we are aware of is [SaplG, pp. 10, Footnote 1], where no details are provided. Therefore, we describe the 
procedure to overcome this for completeness. 

®Note that while A depends on n already, we feel that it is clearer to explicitly designate n as a parameter of K, 
as seen e.g. in Proposition 17. 
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4.2 Sparse interpolation in an interpolation-friendly basis 

In this section we perform deterministic sparse polynomial interpolation in an interpolation-friendly 
basis. The idea is to combine the structure of the Barvinok-Fomin algorithm with some ingredients 
from the Klivans-Spielman algorithm. 

We first briefly review the idea of the Klivans-Spielman algorithm [KSOl, Section 6.3]. Suppose 
we want to interpolate / G Z[xi,... ,Xn] of degree d with m monomials. Their algorithm makes 
use of the map 

(/>c(xi, ... ,Xn) = (2) 

where c = (ci,..., Cn) G (Z+)"'. If c satisfies the property: Ve / e' G /, (c, e) / (c, e'), then we can 
reduce interpolation of multivariate polynomials to interpolation of univariate polynomials: first 
apply the univariate polynomial interpolation (based on the Vandermonde matrix) to get a set of 
coefficients. Then to recover the exponents, modify 0c as 

0(,(xi,...,xn) = (3) 

where pi's are distinct primes, and get another set of coefficients. Comparing the two sets of 
coefficients we can compute the exponents. Note that the components of c need to be small for the 
univariate interpolation to be efficient. To obtain such a c, Klivans and Spielman exhibit a small ~ 
polynomial in n, m, d and an error probability e G (0,1) - set of test vectors c s.t. with probability 
1 — e, a vector from this set satisfies the above property. Furthermore, the components of these 
vectors are bounded by 0{m?nd/e). Their construction was reviewed in Lemma 6, Section 2. 

Now suppose M is a {K, LA)-friendly basis for R C Q[x], and we want to recover / = 
of degree < d, and |F| = m. To apply the above idea to an arbitrary I-friendly basis M, the natural 
strategy is to extract the leading monomials w.r.t. La- However, as each basis polynomial can 
be quite complicated, there are many other non-leading monomials which may interfere with the 
leading ones. Specifically, we need to explain the following: 

(1) Whether extremely large coefficients appear after the map (pc, therefore causing the univariate 
interpolation procedure to be inefficient? 

(2) Whether some leading monomials are preserved after the map 0c? (That is, will the image 
of every leading monomial under 0c be cancelled by non-leading monomials?) 

It is immediate to realize that I-friendly bases are designed to overcome the above issues. (I) 
is easy: by the K-bounded property, for any monomial in /, the absolute value of Coeff(u,/) 
is bounded by K ■ (^^ |oa|)- Therefore, the coefficients of the image of / under 0c are bounded 
by • K ■ I®a|))- (2) is not hard to overcome either; see the proof of Theorem 9 below. 

These properties are used implicitly in the Barvinok-Fomin algorithm. 

There is one final note: if, unlike in the monomial basis case, the procedure cannot produce 
all leading monomials at one shot, we may need to get one t\ and its coefficient, subtract that off, 
and recurse. This requires us to compute tx's efficiently. As this is the property of tx, not directly 
related to the interpolation problem, we assume an oracle which takes an index A G A and an input 
a G N”, and returns txisi)- 

Theorem 9. Let M = {t;, | A G A} 6e a {K, La)- friendly basis for R C Z[x], K = K{\,n), A a 
(0,1) invertible matrix. Given an access to an oracle O = 0(A,a) that computes basis polynomial 
tx{3i) for a G N”, there exists a deterministic algorithm that, given a black box containing f = 
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X^Aer Qa / 0 G Z with the promise that deg(/) < d and |r| < m, computes such an expansion 
of f in time poly {n,d,m,log{Ylx\ax\),log K). 

Proof. We first present the algorithm. Recall the maps (j)c and cj)'^ defined in Equation (2) and 
Equation (3). 

Input: A black box B containing f G R C Z[xi,... ,Xn] with the promises: (1) deg(/) < d] (2) / 
has < m terms in the M-basis. An oracle O = 0(A, a) computing tx{^) for a G N"". 

Output: The expansion / = X^Aer^A^A- 

Algorithm: 1. By Lemma 6, construct the Klivans-Spielman set KS = KS(m,n, 1/3,nd,p). 

2. Eor every vector c in KS, do: 

(a) /c 0. i ■(— 0. d ■(— Ac. 

(b) While i < m, do: 

i. Apply the map (/>d to B — /c (with the help of O), and use the univariate 
interpolation algorithm to obtain g{y) = X]i=o If = d) break. 

ii. Apply the map 4>'^ to B — fc (with the help of O), and use the univariate 
interpolation algorithm to obtain g'{y) = '^\=Qb[y'‘. 

hi. Erom and 6^, compute the corresponding monomial and its coefficient Oe- 
Erom X®, compute the corresponding label u G A, and set ■(— Ug. 

iv. /c ^ /c + CLutu- i^ i + 1. 

3. Take the majority of fc over c and output it. 

We prove the correctness of the above algorithm. As before, let ex be the leading vector of tx 
w.r.t. La. Note that as is (0, l)-matrix, entries in the vector A^ex are in {0,1,... ,nd}. By 
the property of KS(m, n, 1/3, nd,p), no less than 2/3 fraction of the vectors from c G KS satisfy 
that {c,A^ex) are distinct over A G T; call these vectors “distinguishing.” We shall show that for 
any distinguishing vector, the algorithm outputs the correct expansion, so step (3) would succeed. 

Fix a distinguishing vector c. As ex ^ ex' for X ^ X' and A is invertible, there exists a unique 
u G A s.t. e,y achieves the maximum at (c, A^ca) over A G T. As (c, A^oa) = (Ac, ca) = (d, ca), 
by the definition of M being Ly^-compatible, we know that within each tx, d* achieves the unique 
maximum at ca over Et^, the set of all exponent vecors in tx. Thus over all e G /, d* achieves the 
unique maximum at e^. This means that y^^’^A is the monomial in g{y) of maximum degree, and 
could not be affected by other terms. As Coeff(ey, = 1, Coeff(ei,,/) is just the coefficient of t^ 
in the expansion of /. So we have justified that from Step (2.b.i) to (2.b.iii), the algorithm extracts 
the monomial of maximum degree in g{y), computes the corresponding coefficient and exponent in 
/, and interprets as a term axtx- 

To continue computing other terms, we just need to note that c is still distinguishing w.r.t. (/— 
some of the terms within /). This justifies Step (2.b.iv). 

To analyze the running time, the FOR-loop in Step (2) and the WHILE-loop in Step (2.b) 
take 0{mf‘n) and m rounds, respectively. In the univariate polynomial interpolation step, as the 
components in c are bounded by 0{w?‘n ■ nd) and A is (0,1), the components in d are bounded by 
0{m‘^n^d). It follows that k = deg{g) = (d,ei,) = 0{w?‘Ad'^). By the K-bounded property, the 
coefficients of g{y) are of magnitude • A • (X^^,^ |oa|)) • So the running time for the univariate 

interpolation step, and therefore for the whole algorithm, is poly(m, n, d, log(X(A I^ADjlogK). □ 
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5 Schubert polynomials form an interpolation-friendly basis 


In this section, onr ultimate goal is to prove Proposition 17 and Proposition 18, which establish 
that Schubert polynomials form an I-friendly basis. The main theorem Theorem 1 follows immedi¬ 
ately. For this, we need to review some properties of Schur polynomials, the definition of Schubert 
polynomials via divided differences, and the transition formula [LS85]. The transition formula is 
the main technical tool to deduce the properties of Schubert polynomials we shall need for Propo¬ 
sition 17 and Proposition 18. These include Lemma 14 which helps us to find the matrix A needed 
for the Lyi-compatible property, and an alterative proof for Schubert polynomial in ^^P. 

Schur polynomials. For a positive integer £, the complete symmetric polynomial h^(x) G Z[x] is 
the sum over all monomials of degree £, with coefficient 1 for every monomial. We also define ho(x) = 
1, and /i^(x) = 0 for any £ < 0. For a partition a = (oi,..., an) in N”, the Schur polynomial Sq(x) 
in Z[x] can be defined by the Jacobi-Trudi formula as Sq,(x) = det[/iQ,._i+j(x)]jNote that 
deg(sQ(x)) = |a|. Via this determinantal expression, we have 

Proposition 10 ([BF97, Sec. 2.4]). For a G Sa(a) can he computed using 0(|ap • n + n^) 
arithmetic operations, and the bit lengths of intermediate numbers are polynomial in those of a. 

Littlewood’s theorem shows Schur polynomials have positive coefficients. We also need the 
following bound on coefficients - the Kostka numbers - in Sa(x). 

Proposition 11 ([BF97, Sec. 2.2]). For any e G N"", 0 < Coeff(e,Sc(x)) < ^jaj!. 

Definition of Schubert polynomials via divided differences. We follow the approach in 
[LasOS, Las08]. For i G [n — 1], let Xi be the switching operator on Z[xi ,..., Xn\'- 

(xi, . . . , Xj, , . . . , Xn) • — /(^l) • • • ) ^i+l i^i, , Xn). 

Then the divided difference operator di on Z[xi,... ,Xn] is di(/) := 

Definition 12. For v = (ui,... ,n„) G N"", the Schubert polynomial Vv S Z\xi,... ,x„] is defined 
recursively as follows: 

1. If V is dominant, then = x^^x^ ... x'^. 

2. If Vi > Vi+i, then = diY^ where v' = (m,..., Vi+i,Vi -1,... ,Vn). 

It is not hard to see that this defines IV for any v. We list some basic facts about Schubert 
polynomials. 

Fact 13 ([ManOl]). 1. If'v = {vi,... ,Vn) is anti-dominant with k being the last nonzero index, 
Yv equals the Schur polynomial Sa(xi ,... ,Xk) where a = {vk,... ,vi). 

2. As a special case of (1), i/v = (0,..., 0, rc, 0,... 0) where w is at the kth position, then Vv(x) 
is the complete homogeneous symmetric polynomial hw{xi,... ,Xk). 

3. Ifv= [ui,... ,nfc,0,... ,0] G N”, then W, e '^[xi, ...,Xk]. 


12 



The transition formula and its applications. Given a code v G let k be the largest t s.t. 
vt is nonzero, and v' = [ni,..., — 1,0,..., 0]. For convenience let a = (v'). Then the transition 

formula of Lascoux and Schiitzenberger [LS85] is: 

hv — ^ ^ Fu) 

U 

where u G N” satisfies that: (i) (u)(T“^ is a transposition Tik for i < k\ (ii) |u| = |v|. Assuming (i), 
condition (ii) is equivalent to: 

(a) a{i) < (7{k)] (b) for any j s.t. i < j < k, 

either a{j) > cr{k), or a{j) < cr{i). (4) 

Let 'I'v be the set of codes with weight |v| appearing in the transition formula for v, and = 
'I'vU{v'}. Any u G 'Lv is uniquely determined by the transposition Tik, therefore, by some i G [A:—!]. 

The transition formula yields the following simple, yet rarely mentioned^*^ property of Schubert 
polynomials. This is the key to show that the Schubert basis is L^-compatible for some appropriate 
A. For completeness we include a proof here. 

Given v and u in N”, v dominates u reversely, denoted as vi>u, if Vn > Un, Vn+Vn-i > Un+Un-i, 
■ Vn-\ - \-V2>Un-\ -h ^2, 4- \-Vi = Un -\ -h Ui. 

Lemma 14. For u G N”, if is in Y^, then v i> u. Furthermore, Goeff(v, Y^) = 1. 

Proof. We first induct on the weight. When the weight is 1, the claim holds trivially. Assume the 
claim holds for weight < w, and consider v G N” with |v| = u; + 1. We now induct on the reverse 
dominance order, from small to large. The smallest one with weight tc + 1 is [re + 1, 0,... , 0]. As 
F[u;+i,o,...,o] = 2 :)“'''^, the claim holds. 

For the induction step, we make use of the transition formula. Suppose k is the largest i 
s.t. Vi is nonzero, v' = [ui,..., — 1,0,... , 0], and a = (v'). Then by the transition formula, 
Fv = XfcFv' + X^uFu, where u G N"’ satisfies that: (i) is a transposition Tik for i < k\ (ii) 

|u| = |v|. Assuming (i), the condition (ii) is equivalent to that: (a) (T{i) < a{k)-, (b) for any j s.t. 
i < j < k, either a{j) > (j{k), or a{j) < (j{i). Thus u and v' can only differ at positions i and k, 
and Ui > v[ = Vi, Uk Yv'^ < Vk- It follows that v i> u, and it is clear that v 7 ^ u. By the induction 
hypothesis on |v|, each monomial in Y^/ is reverse dominated by v', and Goeff(v', Fv,/) = 1 . As Y^/ 
depends only on xi,..., by 13 (3), each monomial in x^F^' is reverse dominated by v. By the 
induction hypothesis on the reverse dominance order, every monomial in F^ is reverse dominated 
by u, thus is reverse dominated by v and cannot be equal to v. Thus Goeff(v, Fv) = 1, which is 
from XfcFv'. This finishes the induction step. □ 

We then deduce another property of Schubert polynomials from the transition formula. Starting 
with vq, we can form a chain of transitions vq —)■ vi —t- V 2 —)■•••—)• Vj —>•... where Vj G 
The following lemma shows that long enough transitions lead to anti-dominant codes. 

Lemma 15 ([LS85, Lemma 3.11]). Let vq —)■ vi —)••••—)• he a sequence of codes in N"", s.t. 
Vj G i G [^]. If none of Vi’s are anti-dominant, then I < n ■ jvoj. 

^°The only reference we know of is in [Lasl3, pp. 62, Footnote 4]. 
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Based on Lemma 15, we have the following corollary. Recall that for v G N”, is the collection 
of codes (not necessarily of weight |v|) in the transition formula for v. 

Corollary 16. Let vq —)• vi —)■•••—)• V£ be a sequence of codes in N”, s.t. Vj G i G [i]. If 

none ofvi’s are anti-dominant, then i <n - (|vop + |vo|). 

Proof. For w G [|vo|], let i,,, be the last index i in [1] s.t. Vj is of weight w. As none of Vj’s are 

anti-dominant, by Lemma 15, iw-i — iw^n-w-\-l<n - |vo| -|- 1. The result then follows. □ 

In fact, by following the proof of Lemma 15 as in [LS85], it is not hard to show that in Corol¬ 
lary 16 the same bound as in Lemma 15, namely i < n ■ |vo|, holds. 

From Corollary 16, a ffP algorithm for Schubert polynomials can also be derived. 

An alternative proof of Corollary 3. Consider the following Turing machine M: the input to M is 
a code v G N” in unary, a point a = (oi,..., a„) G in binary, and a sequence s of pairs of indices 
{si,ti) G [n] X [n], Si < U, and i G [1] where i := n ■ (|vp -|- |v|). The output of M is a nonnegative 
integer. Given the input (v, a, s), M computes a sequence of codes vq —)■ vi ^ v^, and keeps 

track of a monomial —)■ —)■ • • • —)> x®^, where e* G The pair of indices (sj+i, tj+i) is used 

as the instruction to obtain Vj+i from Vj, and Oj+i from e*. 

To start, vq = v, and e = (0, ...,0). Suppose at step i, e* = (ei,...,en), and Vj = 

(ui,..., Ufc, 0,..., 0), Ufc 7 ^ 0. (A; is the maximal nonzero index in Vj.) 

If Vj is anti-dominant, then equals to some Schur polynomial by 13 (1). Using Proposition 10 
M can compute the evaluation of that Schur polynomial on a efficiently, and then multiply with 
the value OieM output. Note that as will be seen below, the weight of ej is < n • |vop, 

so the bit length of Oiefn] polynomial in the input size. 

In the following Vj is not anti-dominant. M checks whether U+i = k. If not, M outputs 0. 

In the following tj_|_i = k. M then checks whether Sj+i = U+i. 

If Sj+i = tj+i, M goes to step i -|- 1 by setting Vj+i = (ui, ■ ■ ■ ,Vk — 1,0,... ,0), and ej+i = 
(ei,..., efc_i, Cfc -I- 1, Cfc+i,..., Cn). 

If Sj+i < U+i, then M tests whether Sj+i is an index in as follows. It first computes the 
permutation cr ;= ((ui,..., Vk-i,Vk — 1,0,..., 0)) G Sjy, using the procedure described in Section 2. 
Note that N = maxjg[„]{ui -|- i} < |vj| + n < |vo| -|- n. Then it tests whether Sj+i is in using 
Equation (4). If Sj+i is not in then M outputs 0. Otherwise, M goes to step i -|- 1 by setting 
Vj+i — ), and ej_|_i — ej. 

This finishes the description of M. Clearly M runs in time polynomial in the input size. M 
terminates within i steps by Corollary 16. M always outputs a nonnegative integer as Schur 
polynomials are polynomials with positive coefficients. 

Finally note that Fv(a) is equal to ^gM(v, a, s), where s runs over all the sequences of pairs 
of indices as described at the beginning of the proof. By the discussion on in Section 2, this 
puts evaluating on a in ffP. □ 

The Schubert basis is interpolation-friendly. Now we are in the position to prove that the 
Schubert basis is interpolation-friendly. 

Proposition 17. {1^ S Z[xi,... ,Xn] \ v G N”} is K-bounded for K{'v,n) = n^”'(|vp+|v|) . y^jVjT. 
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Proof. The alternative proof of Corollary 3 for Schubert polynomials implies that can be written 
as a sum of at most polynomials /, where / is of the form x® • Sq, \a\ + |e| = |v|. 

The coefficients in Schur polynomial of degree d are bounded by yfd). by Proposition 11. The claim 
then follows. □ 

Proposition 18. {Tv £ ■ ■ ■ ,Xn] | v G N”'} is LA-compatible for A = Xij = 0 if 

i < j, and 1 otherwise. The leading monomial ofY^ w.r.t. La is x''. 

The matrix A in Proposition 18 is the lower triangular matrix of I’s on the diagonal and below. 
Compare with that for Schur polynomials, described in Section 4.1. 

Proof. This follows easily from Lemma 14: note that for any c = (ci,... ,Cn) G u G Tv, 

(^C,u) = Ci(ui H- Un)-\-C 2 {U 2 H-h Un) H-h Cn^n < Ci(ui H-hUn) + 02(^2 H-+ 

• • • + CnVn = (^c, v). As Cj > 0, the equality holds if and only if u = v. □ 

Now we conclude the article by proving the main Theorem 1. 

Proof of Theorem 1. Note that -y^jvIT is upper bounded by 

and recall |v| = deg(Tv). Then combine Proposition 17, Proposition 18, and Theorem 9. □ 
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A An alternative proof of skew Schubert polynomials in VNP. 

By Equation Equation (1), skew Schubert polynomial can be written as: 

Pr/cr(^) ^2; • • • ^ n ') 

N-1 

:= Y. = E n (5) 

C C i=l 

Note that ew = 0. 

Let m be the length of an increasing chain. The first and second indices of each label in an 
increasing chain are encoded by V x m 0 — 1 matrices g and b respectively. In these matrices the 
variables are listed along the row and the edges of a chain along the column. So for each column, 
there will be a 1 in the row which corresponds to the index that appears in that label. Thus in 
case of the g matrix it indicates which variable should be multiplied in the monomial. 

All permutations are represented hy N x N 0 — 1 matrices, acting on length-V column vectors. 
Ho and V, representing permutations a and vr respectively are given. Let Wi,.. .Wm be the 
intermediate permutations in the increasing chain. 

Now consider the following polynomials. 

N m N 

hiN = HE ^k^dkj) (6) 

i=l j=l k=l 

hiN encodes the monomial for a given increasing chain, which depends on g. 

In the following we shall gradually build up a series of polynomials, which basically characterize 
the property that g and b form an increasing chain. 


m N N 

h2N =hk n (1 - mij(wt)ik)) • ( 7 ) 

t=l i,j,l,k i=l j=l 

where the second product is over all I < i,j,l,k < N such that i = I iS j ^ k. h 2 N encodes valid 
permutation matrices. That is, it is non-zero iff Wi,... Wm are valid permutation matrices, that 
is each row and column contains exactly one 1. 


h^N [(1 9ik9kj)] ■ [(1 


( 8 ) 


where both the products are over all 1 < /c < m and 1 < j < such that i ^ j. is non-zero 
iff each column of g and b has at most one 1. 

m N m N 

^4Ar = [H X] ■ [11 

j=l k=l j=l k=l 

h^N is non-zero iff there is at least one 1 in each column of g and b. 

N 

h5N = n [1 - (i^rnh " Vij)] ■ [1 + - Vi,)] (10) 

*J = 1 

/i 5 Ar is non-zero iff Wm = V. 

N 

hm = n [1 - - (.Wt-lTij)ab)Ul + {{Wt)ab - {Wt-lTij)ab)] ' bjt (h) 

a,6=1 

h^ijt is non-zero iff for a particular transposition {i,j) and for a pair of consecutive permutations 
VFf_i and Wt-, (a) Wt = Wt-iTij and (b) second index of the label is given according to definition 
of labeled Bruhat order. 

N N 

hrt = EE (l^t—l)ifc ■ (ffh—l)jZ ■ [1 (bhf—l)af)] ■ Agjji (12) 

ij=l k,l=l i<a<j 

k<l k<b<l 

h^t is non-zero iff for any pair of consecutive permutations a', vr' being encoded by Wt-i and Wt 
respectively such that vr' = r^cr' , we have cr'{i) < cr'{j) {i < j) and for every i < k < j, either 
a'{k) < a'{i) or cr'{k) > cr'{j). That is, a' < vr' in the Bruhat order. 


hsN — ^ gst ( 13 ) 

t=l 

where 1 < s < A^. h^N is non-zero iff the sequence of permutations is correct, namely, maintaining 
the Bruhat order so that the labels are given accordingly. 

m-l N N N N 

hgN = BE 9it[ ^ ^ 9j,t+i + 9i,t+i ( ^~~^ bkt ■ ^ ^ ^z,t-i-i)] (1^) 

t=l i=l k=l l=k 

hgjy is non-zero iff each pair of labels in a chain respects the increasing lexicographic order. 

We define the polynomial : 

/iAr(xi, ...XN,g,b,Wi,... Wm) 

= hiN ■ h2N ■ hsN ■ hitN ■ h^N ■ hsN ■ hgN ( 15 ) 
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For each assignment to g, b, Wj, is either 0, or a monomial corresponding to one correct increasing 
chain. It is clear that the size of a straight line program to evaluate is polynomial in N. So 
is in VP. 

The skew Schubert polynomial, then can be given by 

.. .xn) = ^ hNixi,...XN,g,b,Wi,...Wm) (16) 

g,b,Wi,...Wm 

g and b have mN elements and since m is 0{N‘^) so each has 0{N^) elements. Each Wi has 
entries. Thus the summation is over 0-1 strings of polynomial length. 

This proves is p-definable and hence is in VNP. 
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